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ABSTRACT 

We explore the dependence of pre-supernova evolution and supernova nucleosynthesis yields on the 
uncertainties in helium burning reaction rates. Using the revised solar abundances of Lodders (2003) 
for the initial stellar composition, instead of those of Anders k. Grevesse (1989), changes the supernova 
yields and limits the constraints that those yields place on the ^^C(a,7)^^0 reaction rate. The 
production factors of medium-weight elements (A = 16-40) were found to be in reasonable agreement 
with observed solar ratios within the current experimental uncertainties in the triple-a reaction rate. 
Simultaneous variations by the same amount in both reaction rates or in either of them separately, 
however, can induce significant changes in the central ^^C abundance at core carbon ignition and in the 
mass of the supernova remnant. It therefore remains important to have experimental determinations 
of the helium burning rates so that their ratio and absolute values are known with an accuracy of 10% 
or better. 

Subject headings: Nucleosynthesis, Nuclear Reactions, Sun: Abundances, Supernovae: General 



1. INTRODUCTION 

Throughout the past three decades much experimental 
and theoretical effort has been dedicated to determining 
the rate of the ^^C{a, j^^O reaction. It, and the triple-a 
(3a) reaction, are responsible for both energy generation 
and the nucleosynthesis of C and O during stellar helium 
burning; the ratio of their rates determines the ratio of 
carbon to oxygen at the completion of core helium burn- 
ing. This ratio, in turn, strongly influences the subse- 
quent evolution of Type II supernova (SNII) progenitors 
(> 9 Mq stars) , affecting both the pre-supernova stellar 
structure and the post-explosive nucleosynthesis. 

Although progress has been achieved in the labora- 
tory evaluation of the ^^C(a,7)^^0 reaction rate, Ra,i2, 
there arc significant uncertainties in its extrapolation to 
the low energies relevant to hydrostatic helium-burning 
in stars (~ 300 keV). In a recent review, Buchmann & 
Barnes (2006) recommend S(300 keV) = 145 kcV b with 
errors in the range of 25%-35% for the measured astro- 
physical S-factor of ^'^C{a,j)^^0. Stellar evolution cal- 
culations have shown (Weaver & Wooslcy 1993; Boyes 
et al. 2002 [reported in Woosley et al. 2003; Woosley 
& Heger 2007]) that such an uncertainty has major ef- 
fects on SNII nucleosynthesis and on the mass of the 
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pre-collapse core. 

These calculations examined the changes in the pro- 
duction factors (defined as the ratio of the average iso- 
topic mass fraction of nuclides in the ejecta to their solar 
mass fraction) induced by varying Ra.i2- For these stud- 
ies, however, the pre-supernova isotopic mass fractions 
were used in determining the production factors. Under 
the crude assumption that SNII progenitors of close-to- 
solar metallicity are the main contributors to the ob- 
served solar abundances for the medium-weight isotopes 
(A — 16-40), very similar production factors arc desirable 
for those nuclides. The reaction rate that produces the 
smallest spread in production factors was found; Boyes 
et al. (2002) find the smallest spread for a narrow range 
in ^^C(a,7)^®0 reaction rates that is only 10% wide. 

The calculations assumed a fixed value of the 3a re- 
action rate, Rsa- This is a reasonable assumption since 
i?3Q has significantly smaller experimental uncertainties, 
about 10%-12% (Tur et al. 2006; Austin 2005). However, 
if an accuracy of 10% in the ratio Rsa/ Ra,i2 is required, 
the present accuracy of the 3a rate is insufficient. 

Besides uncertainties in nuclear reaction rates, uncer- 
tainties in the initial isotopic composition of stars affect 
their evolution and nucleosynthesis. Most recent stud- 
ies of SNII evolution used the abundances of Anders & 
Grevesse (1989); to our knowledge there are no sys- 
tematic studies using the more recent abundance set of 
Lodders (2003). 
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Fig. 1. — Lines and dots show the three sets of simulations done 
for eaeh star of a given initial mass and a given initial solar abun- 
dance distribution. For the blue squares, R^ci is held constant (at 
its value by Caughlan & Fowler 1988), and i?c<,i2 is varied (A). For 
the red dots, both reaction rates are varied by the same percentage, 
so their ratio remains constant (B). For the green crosses, Rc<,i2 
is held constant at 1.2 times the rate recommended by Buchmann 
(1996) and R^ci is varied (C). 



In this paper we describe an extensive set of calcula- 
tions to determine how SNII nucleosynthesis and other 
stellar properties vary when R^a and Ra,i2 arc varied. 
These calculations were repeated for the two abundance 
sets: Anders & Grevesse (1989) and Lodders (2003). We 
repeated the calculations of Boyes et al. (2002) to ensure 
that any small changes in procedures are unimportant. 
Another improvement is that the results of Boyes et al. 
(2002), as weh as those of Weaver & Woosley (1993), 
were based on pre-SN nucleosynthesis, but some of the 
abundances examined are known (Weaver & Woosley 
1993; Woosley et al. 2002) to be modified in the SN ex- 
plosion. The simulations presented in this paper include 
explosive nucleosynthesis. 

A description of the stellar models and the range of the 
calculations is given in Section 2. The differences in the 
stellar structure and nucleosynthesis resulting from dif- 
ferences in solar abundance sets are presented in Section 
3. In Section 4, we compare the stellar evolution impli- 
cations of the uncertainties in the 3a and ^^C(a,7)^^0 
rates. 

2. COMPUTED MODELS 

Stars with initial masses from 13 to 27 were con- 
sidered. All models were calculated using the implicit, 
one dimensional, hydrodynamical stellar evolution code 
KEPLER. Since its first implementation in 1978 (Weaver 
et al. 1978) the code has undergone several major revi- 
sions with improvements to the physical modeling of the 
stellar structure and to the nuclear reaction networks 
(Woosley & Weaver 1995; Rauscher et al. 2002; Woosley 
ct al. 2002). A small network directly coupled to the stel- 
lar model calculation provides the approximate nuclear 
energy generation rate. A larger "adaptive" network is 
used to track nucleosynthesis. The large network auto- 
matically adjusts its size to accommodate the current nu- 
clear flows and progressively grows from several hundred 
isotopes during hydrogen burning to more than 2000 iso- 



topes at explosive burning. The treatment of convection, 
semi-convection, and overshoot mixing is as described in 
Woosley & Weaver (1988) and Woosley et al. (2002). We 
do not include the effects of rotation and magnetic fields. 
Stars are first evolved from the zero age main sequence 
to pre-supernova, i.e., from central hydrogen burning to 
iron core collapse, and are then exploded. The explosion 
is parameterized by a piston at a constant Lagrangian 
mass coordinate and has two important specifications: 
its location in mass (i.e., the initial mass cut) and the 
total kinetic energy of the ejecta at infinity (here, 1 yr 
after the explosion). See Woosley & Heger (2007) for a 
more complete description of the piston parameters. 

The values of these two parameters are chosen to fit 
reasonably well within the range of observational con- 
straints. The explosion energy was set to 1.2B (B for 
Bethe, IB = 10^^ ergs). Supernova 1987A is thought 
to be an 18-20 Mq star which exploded with an esti- 
mated energy in the range 0.6-1.5B with an uncertainty 
of perhaps 50% based on the observed light curve and 
velocity (Arnett et al. 1989). The initial mass cut was 
placed at the base of the oxygen burning shell, a loca- 
tion associated with a large density drop, and hence dy- 
namically important to generating successful explosions 
(Janka 2007). Specifically, we chose to place the piston 
at the location in the star where the entropy S reaches a 
value of S* = 4fcB/baryon (Woosley & Heger 2007), be- 
yond which a large rise in entropy, hence a drop in den- 
sity, is observed. The piston location cannot be below the 
surface of the iron core or neutron-rich species in the iron 
group will be overproduced; it cannot be above the base 
of the oxygen shell or typical neutron star masses will be 
too large (Woosley & Heger 2007). Our nucleosynthesis 
studies do take into account all strong and weak reac- 
tions during oxygen shell burning, including the slight 
neutron excess resulting in this burning phase. KEPLER 
calculations by Woosley & Heger (2007) showed that ex- 
plosion energies of either 1.2B or 2.4B and mass cuts at 
the base of the oxygen burning shell or at the edge of 
the iron core gave very similar nucleosynthesis, except 
for the iron peak nuclei. We note, however, that in a 
recent study by Young & Fryer (2006), both elemental 
and isotopic yields beyond silicon were found to be very 
sensitive to the explosion energy. 

Three separate studies were done for stars of 15, 20, 
and 25 for both solar abundance sets, Anders & 
Grevesse (1989) and Lodders (2003): (A) i^g^ was kept 
constant (at its value from Caughlan & Fowler 1988) and 
Ra,i2 was varied; (B) both rates were varied by the same 
factor, so their ratio remained constant; and (C) Ra,i2 
was held constant at 1.2 times the rate recommended by 
Buchmann (1996) and R^a was varied. The ranges of 
those variations are shown in Figure 1. For the Anders & 
Grevesse (1989) abundances, we additionally computed 
stars of 13, 17, 19, 21, 23 and 27 Mq to have a better 
sampling of the initial mass function (IMF) (see eq. 1), 
in order to better integrate over intrinsic star to star vari- 
ations and thereby reduce the impact of numerical noise 
in the production factors. The isotopic mass fractions 
from all the stars in a given study were then averaged 
over an IMF with a slope of 7 = —2.6 (Scalo 1986) and 
divided by their solar mass fraction, giving the produc- 
tion factor of each isotope. The slope 7 is defined by the 
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Fig. 2. — Ratio of the Loddcrs (2003) abundances to the Anders 
& Grevesse (1989) abundances as a function of mass number (up 
to strontium). Isotopes of each element have the same color and 
are connected by lines. 



Table 1. Simulation series 



Label Description 



AAX^ Anders & Grevesse (1989); Ra,i2 varied 
ABX'' Anders & Grevesse (1989); iJs^, Ra,i2 varied 
ACX^ Anders & Grevesse (1989); Rsa varied 

LAX'' Loddcrs (2003); Ra.i2 varied 
LBX'' Loddcrs (2003); Rsq, Ra,i2 varied 
LCX'' Loddcrs (2003); Rsa varied 

''X=2, if IMF average over 2 stars (15 and 25 Mq); X=8, if IMF 
average over 8 stars (13, 15, 17, 19, 21, 23, 25, and 27 Mq) 



equation: 



(1) 



where ^{logM) is the IMF defined in units of the number 
of stars per (base 10) logarithmic mass interval M per 
square parsec of the Galactic disk, M is the initial mass 
of the star in solar masses, and A and 7 are constants 
(Weaver & Woosley 1993). 

We will adopt a three-character notation to label our 
plots, e.g., LA2 (see Table 1). The first character can be 
an L (to denote the Lodders 2003 initial abundances) 
or an A (for the Anders & Grevesse 1989 initial abun- 
dances). The second character denotes the study: A 
when R^a was kept constant and i?a.i2 was varied; B 
when both rates were varied by the same factor, so their 
ratio remained constant; and C when Ra,i2 was held con- 
stant and i?3Q, was varied. The third character is a num- 
ber; it is 2 when the production factors are averaged over 
two stars (15 Mq and 25 Mq) and 8 when the average is 
over eight stars (13, 15, 17, 19, 21, 23, 25, and 27 Mq). 
When no third character is present, no average has been 
performed, as in the case where the numbers only apply 
to a single star. 

3. SENSITIVITY TO DIFFERENCES IN SOLAR 
ABUNDANCES AND REACTION RATES 

The differences in the two recent solar abundance de- 
termination are shown in Figure 2. For Lodders (2003) 



compared to Anders & Grevesse (1989), the abundances 
are the following: for CNO they are lower by about 30%; 
for CI, Kr, Xe, and Hg they are higher by more than 40%; 
and for most other metals they arc higher by about 15%. 
As a consequence, the overall solar mass fractions change 
from Xa = 0.7057, = 0.2752 and Zq = 0.0191 for the 
old set to 0.7110, yo = 0.2741, and Zq = 0.0149 

for the new set. 

3.1. The effect on the production factors 

The study with R^a constant and Ra,i2 varied, is 
an elaboration of two previous studies using the An- 
ders & Grevesse (1989) abundances (Weaver & Woosley 
1993; Boyes et al. 2002). Relative to Weaver & Woosley 
(1993), our models also include mass loss due to stel- 
lar winds, as described in Woosley & Heger (2007). As 
noted above, explosive nucleosynthesis is also included. 
The same study done with the Lodders (2003) solar 
abundances is entirely new and demonstrates the uncer- 
tainties in determining R^a and Ra.i2 using astrophysical 
models. 

Based on SNII nucleosynthesis considerations. Weaver 
& Woosley (1993) predicted an S-factor at 300 keV of 
~170 keV b; or more precisely a rate of 1.7 ± 0.5 times 
that of Caughlan & Fowler (1988). This constrained 
i?a,i2 to a range of about 30%. The same study was re- 
peated later by Boyes et al. (2002) (reported in Woosley 
et al. 2003; Woosley & Heger 2007) with improved stel- 
lar models [newer opacities, added mass loss, finer stellar 
zoning, and finer grid of ^^C{a,j)^^0 rates] and found 
a best fit of 175 keV b or about 1.2 times the value of 
S(300 keV) suggested by Buchmann (1996) (146 keV b). 
This study concluded that Ra 12 needed to be known to 
< 10% (Woosley et al. 2003; Woosley & Heger 2007). 

In Figures 3A, 3C and 3D, we illustrate our results 
for the Anders & Grevesse (1989) abundances. Figure 
3A shows the production factors averaged over two stars 
(AA2), and their rms deviations for the same set of iso- 
topes as those selected by Boyes et al. (2002). Figure 
3C does the same for a larger group of stars (AA8) and 
Figure 3D does the same for a larger set of medium- 
weight isotopes (now including ^^F, '^^P, ^^Cl, and ^^^K). 
If SNII are indeed the major site of production of all 
medium- weight elements (A = 16-40), then those ele- 
ments should have similar production factors at a point 
where their rms deviations are minimum. For the Anders 
& Grevesse (1989) abundances, the conclusion seems ro- 
bust; the position of the minimum is well defined at a 
rate of 1.2 times the Buchmann (1996) rate for different 
sets of stars and nuclides, although the details of the rms 
curves vary somewhat. This conclusion agrees with the 
earlier work by Weaver & Woosley (1993) and Boyes et 
al. (2002). 

For the Lodders (2003) initial abundances the results 
are less definitive, as shown in Figure 3B. The average 
production factors at the minimum arc about the same 
for both abundance sets. However, the rms curve now 
has a much broader minimum, again centered around 1.2 
times the Buchmann (1996) rate, but extending from a 
rate multiplier of 0.9 to 1.5. The spread in production 
factors at the minimum is larger by about a factor of 
2. These production factors apparently provide a much 
less stringent constraint on i?Q,i2, allowing a range of 
±25% around the central value of 1.2 times the Buch- 
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Fig. 3. — Panel A: Production factors for AA2 and their rms deviations from the mean for the same set of isotopes as those selected 
by Boyes et al. (2002). A multiplier of 1 means a rate of 1 times the rate recommended by Buchmann (1996). Panel B: Same as Panel 
A, but for LA2. Panel C: Same as Panel B, but for AA8 Panel D: Same as Panel C, but with the addition of the production factors for 
19p^ 31p^ 35^1^ jjj^j 39^;. 



mann (1996) rate. This is unfortunate, since it means 
that one cannot so strongly limit the uncertainty in Ra,i2 
using SNII calculations of production factors. 

Our results are for post-explosion values of the produc- 
tion factors, whereas the two previous studies stopped at 
the pre-supernova stage. We found that nuclides beyond 
^®Si, such as ^°Ca, ^^Ar and ^^S, were significantly mod- 
ified during explosion, often by a factor of 1.5 or more for 
the Anders & Grevesse (1989) abundances. Yet, these 
modifications did not greatly change the earlier results 
for production factors. 

The 20 M0 star showed a peculiar behavior: large over- 
productions were found for '^^P, ^^Cl, and '^^K for some 
choices of R^a and Ra.i2- This was also observed and ex- 
plained in Rauscher et al. (2002). The over-productions 
are attributed to the merging of the convective oxygen-, 
neon- and carbon-burning shells about 1 day before the 
explosion, thereby carrying neutron sources such as ^^Ne 
and ^^Mg to depths where they burn rapidly and pro- 
vide neutrons for capture reactions. Because of these 
peculiarities, we excluded the 20 M© star from our re- 
sults. 



3.2. Variations in the carbon mass fraction at central 
carbon ignition and in the remnant masses 

We also explored the change of the central carbon mass 
fraction at core carbon ignition and of the remnant mass 
after explosion, both as a function of the initial solar 
abundances of the stars and as a function of variations 
in the helium-burning reaction rates. We illustrate the 
results in Figures 4 (carbon mass fraction), and 5 (rem- 
nant masses). The remnant masses are the gravitational 
masses of the resulting neutron stars or black holes. They 
are based on the baryonic mass below the piston (i.e., the 
mass enclosed within a radius reaching out to the base of 
the oxygen shell at the pre-supernova stage) corrected for 
the binding energy (Zhang et al. 2007) according to the 
approximation given by Lattimcr & Prakash (2001). In 
our study, none of the stars had any significant fallback 
after explosion. 

The variations in the central carbon mass fractions 
are smooth, but we see a very sensitive dependence of 
the remnant masses on the solar abundance set used for 
the initial stellar composition. To disentangle and assess 
the magnitude of the effects compared to observational 
data, however, would require a detailed population syn- 
thesis study of remnant masses as a function of metal- 
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Fig. 4. — Panel A: Carbon mass fraction at the center of the star at core carbon ignition for 15, 20, and 25 stars and for the AA 
series. Panel B: Same as Panel A, but for LA. Panel C: Same as Panel A, but for LC. Panel D: Same as Panel A, but for LB. 



licity which is beyond the scope of this paper. The pre- 
dicted remnant masses are also strongly dependent on 
the precise reaction rates used in the pre-supcrnova evo- 
lution, often varying by 0.2 Mq or more, over ranges of 
±2ct experimental errors of the reaction rates. 

The remnant mass may determine the relative popula- 
tion of neutron stars and black holes resulting from SNII 
explosions. Lattimer & Prakash (2007) have surveyed 
the available data on neutron star masses. Their conclu- 
sions are the following: (1) While some masses in excess 
of 2 Mq have been reported, "it is furthermore the case 
that the 2a errors for all but two systems extend into the 
range below 1.45 Mq, so caution should be exercised be- 
fore concluding that firm evidence of large neutron star 
masses exists." And (2) the smallest "reliably estimated 
neutron star mass is 1.18 ±0.02 Mq" . While these uncer- 
tainties make it difficult or impossible to use neutron star 
masses to place limits on reaction rates or abundances, 
for orientation we have placed lines in Fig. 5, at values of 
1.7 and 2.0 Mq as possible maximum masses for neutron 
stars. 

4. COMPARING CHANGES IN TRIPLE-a AND ^'^C{a,-y)^^0 

RATES 

In this section, we discuss the relative importance of 
the uncertainties in R^a compared to the 2 times larger 



uncertainties in Ra,i2- Figures 6A and 6B show the 
production factors of some medium- weight isotopes (the 
same set as Boyes et al. 2002) as a function of the triple-a 
reaction rate in two of our studies: R^a varied and Ra,i2 
constant, and both reaction rates varied by the same fac- 
tor. The variations in the production factors (Figure 6A) 
over a range of one standard deviation a {3a multiplier 
from 0.88 to 1.12) are small, although there are larger 
deviations for 2a differences. 

We find a very sensitive dependence of the remnant 
masses on the helium burning reaction rates, and on 
the initial solar abundance set used. The smooth de- 
crease in the carbon mass fraction as a function of in- 
creasing i?Q,i2, or decreasing R^a, is expected. The fol- 
lowing argument is commonly given to explain the gen- 
eral increasing trend of the remnant masses when the 
Ra,i2 is increased (seen in Figure 5 A for instance): a 
smaller rate gives a larger carbon abundance after he- 
lium burning. During carbon shell burning, this larger 
abundance supports longer and more energetic burning 
which allows the central regions to cool to lower entropy. 
The lower entropy, in general, gives smaller iron cores 
(hence remnants) for stars of a given main-sequence mass 
(Woosley et al. 2003). Figure 4 also shows that smaller 
stars make more carbon than larger ones, reflecting their 
higher density, and tend to have smaller remnants fol- 
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Fig. 5. — Gravitational mass of the remnant (neutron star or black hole) after explosion for 15, 20, and 25 Mq stars. The dotted lines 
at 1.7 and 2.0 Mq mark possible maximum masses for neutron stars (Lattimer & Prakash 2007); for heavier masses, black holes may be 
formed. Panel A: AA series. Panel B: LA series. Panel C: LC series. Panel D: LB series. 
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Fig. 6. — Panel A: Production factors and their rms deviations from the mean for some medium-weight isotopes (the same set as Boyes 
et al. 2002) for LC2. Panel B: Same as Panel A, but for LB2. 



lowing explosion (as seen in Figure 5), which supports 
the previous argument. When looking at the remnant 
masses for the 25 Mq star (Figure 5C) the same argu- 
ment seems to break down at least partly. One expects a 
general decreasing trend in the remnant masses for higher 



triple-a rates, whereas one sees an increase for a multi- 
plier larger than one. The non-monotonic behavior of 
remnant masses can be understood as a result of the in- 
teraction of subsequent burning shells. This causes the 
behavior of the remnants for the 25 Mq star of Figure 
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5C. 

In Figure 5, variations within the current experimen- 
tal range of uncertainties {2a) of both R^a and Ra,i2 
cause significant changes in the remnant mass. The rem- 
nant mass curves look smooth for the 25 Mq star, but, 
they show an oscillatory behavior with rapid variations 
(over a small rate multiplier range) for our 15 and 20 Mq 
stars. In particular. Figure 5D shows very strong fluctu- 
ations in remnant masses, when the ratio of the helium 
burning reactions is kept constant (LB), despite the very 
smooth change of the carbon mass fractions. These os- 
cillations are likely due to small numerical noise in the 
models originating from temporal and spatial discretiza- 
tion, combined with a sharp transition in the stellar evo- 
lution past helium burning as a function of the carbon 
mass fraction, where an additional burning shell ignites 
or does not ignite beyond a certain threshold. 

These observations lend support to the idea that varia- 
tions in both R^a and Ra,i2 are important, not just their 
ratio or their relative variations. An increase of 10% in 
Rsa gives the same amount of increase in the central car- 
bon mass fraction as an 8% decrease in Ra,i2, in close 
agreement with the findings of Woosley & Heger (2007) 
for a simple calculation at given temperatures and den- 
sities. A 27% decrease in both reaction rates is required 
to produce the same amount of increase in the central 
carbon mass fraction when the two rates are multiplied 
by the same factor. 

5. CONCLUSION 

Our simulations show that multiple uncertainties sig- 
nificantly influence the evolution and nucleosynthesis of 
SNII in current one dimensional massive star and su- 
pernova models. The notable effect of differences in so- 
lar abundance sets is one example. Using the Loddcrs 
(2003) abundances rather than the previous standard 
set by Anders & Grevesse (1989), appears to reduce 
the precision with which SNII simulations of production 
factors can be used to constrain Ra,i2 to ±25%. The pro- 
duction factors of medium-weight elements (A = 16-40) 
were found to be about constant within the current la 
experimental uncertainties in the triple-a reaction rate. 
However, variations within the 2a experimental errors in 
either helium-burning reaction rate do induce strong rms 
deviations for the production factors far from the central 
values of those rates. 

We want to issue a caution, however, about our 
very approximate treatment of galactochemical evolu- 
tion. Stars from different initial metallicities contribute 
to the solar abundance pattern. Here we took the ap- 
proximation that the stars which contributed most are 
those of about solar initial abundance, within roughly a 
factor of 2. Although we did not try to obtain a pre- 
cise quantification of the uncertainties due to the form of 
the initial mass function (IMF), the results of our study 
were not changed in any significant way by substitut- 
ing a Salpeter IMF for the Scalo IMF used throughout 
this study. Another physics uncertainty which could af- 
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feet the pre-supernova structure and supernova nucle- 
osynthesis yields is the treatment of hydrodynamics in- 
cluding convection and boundary layer mixing such as 
overshoot and semi-convection. These uncertainties have 
been shown (Woosley & Weaver 1988; Young et al. 
2005) to have effects comparable to uncertainties in nu- 
clear reaction rates, for instance, regarding predictions 
of both carbon mass fraction and remnant mass. One 
more issue concerns the poorly understood interactions of 
burning shells. These effects were discussed in Rauscher 
et al. (2002), and we have pointed out above how they can 
affect nucleosynthesis for a 20 Mq star. Such effects have 
also been confirmed in multi-dimensional calculations of 
pre-supernova stars (Meakin & Arnett 2006). The effects 
of uncertainties in the calculation of mass loss and the 
possible effects of a binary companion could also be im- 
portant. It would be useful to have a numerical estimate 
of the implications of all these uncertainties. However, 
to make something better than a guess would involve a 
suite of calculations much larger than the already exten- 
sive set we have performed. Eventually, perhaps, these 
effects will be sufficiently well known to permit a reliable 
estimate of overall uncertainties. However, even then the 
large effects of uncertainties in the nuclear reaction rates 
will likely remain. 

Within the scope of our study, uncertainties within the 
current errors in the rates of the helium burning reac- 
tions, both triple-a and ^^C{a,j)^^0 have been found 
to induce strong changes in the remnant mass of mas- 
sive stars, highlighting the fact that those rates are inde- 
pendently important. The changes in remnant mass can 
have consequences for the typical neutron star masses. 
Hence, determining the helium-burning reaction rates is 
an essential ingredient to the theoretical understanding 
of the populations of neutron stars and black holes. 

Taken together, our results for SNII evolution support 
the need for improved measurements of both the helium- 
burning reaction rates, with the goal that their ratio is 
known to within 10%. This is particularly important if 
predictions of average remnant masses are to be reliable. 
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